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1.  Introduction 


In  recent  investigations  [2],  [3],  [9],  Banks,  Inman,  and  their  colleagues  carried  out  a  series 
of  experimental  and  computational  investigations  on  linear  damping  models  for  composite 
material  flexible  structures.  Other  investigations  [1],  [7]  involving  experiments  with  flexible 
structures  point  to  the  need  to  understand  nonlinear  aspects  (in  particular,  damping)  of 
flexible  structures.  In  this  paper,  we  present  a  first  step  toward  development  of  a  rigorous 
theoretical  foundation  for  computational  methods  to  study  nonlinear  mechanisms  in  Buch 
structures. 

Traditionally,  nonlinearities  in  structural  models  usually  arise  when  one  needs  to  con¬ 
sider  large  amplitude  displacements.  However  a  rather  substantial  engineering  literature 
(see  Chap.  7  of  [19])  exists  on  a  wide  variety  of  “nonlinear”  effects  in  structures.  There 
are  numerous  sources  proposed  for  these  nonlinearities;  they  can  be  geometric  (nonlinear 
stretching;  large  curvature)  or  material  (nonlinear  stress-strain  laws;  nonlinear  damping) 
in  nature.  In  mathematical  models  these  can  be  manifested  in  nonlinear  boundary  con¬ 
ditions,  nonlinear  stiffness  operators,  and/or  nonlinear  damping  operators.  For  example, 
in  the  usual  derivations  [11],  [20],  [24]  for  the  Euler- Bernoulli  theory  of  beams  or  the 
Love-Kirchoff  theory  of  plates,  one  must  make  several  “linearizing”  assumptions.  The 
commonly  postulated  Hooke’s  law  is  a  linear  stress-strain  constitutive  law.  Even  under 
this  assumption,  one  has  a  nonlinear  moment- curvature  relationship  (and  hence  a  non¬ 
linear  stiffness  operator  in  the  resulting  partial  differential  equation  model)  unless  one 
makes  the  usual  approximations  (e.g.,  see  Chap.  11  of  [20])  in  the  curvature  expression 

<l  fry 

l/R  =  uxxj  [l  +  u2x ]  arising  in  elementary  beam  theories. 

While  the  ideas  presented  below  can  be  used  to  develop  a  methodology  to  study  nonlinear 
stiffness  operators  (which,  of  course,  can  be  of  great  importance  in  motions  involving  large 
geometries  or  nonlinear  material  properties),  we  shall  focus  here  on  a  theory  for  nonlinear 
damping  operators.  Our  approach  will  allow  one  to  study  damping  mechanisms  that  do  not 
decouple  across  the  modes  of  the  structure  and  this  has  substantial  practical  significance. 
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The  efforts  of  Balakrishnan,  Taylor,  and  their  colleagues  (see  [1]  and  the  discussions  and 
references  therein)  have  produced  evidence  that  in  complex  structures  such  as  the  SCOLE 
structure  at  NASA  Langley,  one  cannot  use  linear  damping  models  where  the  damping 
decouples  across  the  “linear  undamped  modes”  to  successfully  describe  experimental  ob¬ 
servations,  The  presentation  in  [1]  suggests  that  modal-based  nonlinear  damping  might 
be  adequate.  However,  in  light  of  the  results  in  [3],  [9],  where  it  is  found  that  for  even 
simple  structures  one  can  be  in  the  unfortunate  situation  where  damping  does  not  decou¬ 
ple  across  the  modes,  this  issue  is  far  from  settled.  Indeed,  it  is  not  clear  whether  such 
damping  should  be  modeled  with  a  nonlinearity  or  whether  it  might  be  adequately  mod¬ 
eled  with  a  linear  mechanism  that  does  not  decouple  (i.e,  cannot  be  described  via  modal 
damping  ratios  in  the  usual  engineering  manner).  Whatever  the  case,  it  is  clear  that  fur¬ 
ther  investigations  are  necessary.  Thus  one  requires  methods  to  study  linear  and  nonlinear 
damping  mechanisms  in  models  of  second  order  systems  where  one  does  not  assume  that 
the  damping  model  can  be  simplified  by  modal  representations. 

Our  goal  then  is  to  develop  such  a  methodology  for  nonlinear  systems.  The  theoret¬ 
ical  basis  for  the  linear  studies  in  [2],  [3],  [9]  was  an  abstract  approximation  framework 
developed  by  Banks  and  Ito  in  [4],  It  is  our  purpose  here  to  develop  a  theory  for  nonlin¬ 
ear  second  order  systems  that  parallels  that  of  [4].  Since  we  shall  use  monotone  operator 
theory  and  nonlinear  evolution  systems  (a  theory  which  itself  has  been  shown  to  possess 
substantial  limitations  in  treating  very  general  nonlinearities),  we  expect  this  to  be  only 
an  initial  contribution  from  our  studies,  However,  the  theory  presented  below  will  allow 
us  to  treat  a  large  class  of  nonlinear  damping  operators  that  are  thought  to  be  of  practical 
importance. 

We  develop  an  abstract  approximation  framework  and  convergence  theory  for  the  identi¬ 
fication  of  nonlinear  dissipation  or  damping  mechanisms  in  second  order  distributed  param¬ 
eter  systems.  More  precisely,  we  consider  the  estimation  of  nonlinear  maximal  monotone 
damping  operators  in  abstract  wave  equations  with  linear  symmetric  stiffness  operators. 
Our  treatment  here  is  in  the  spirit  of  our  earlier  work  on  first  order  systems  in  [8]  and 
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is  an  extension  of  Banks’  and  Ito’s  results  in  [4]  concerning  the  identification  of  second 
order  linear  systems,  By  rewriting  the  underlying  second  order  initial  value  problem  as 
an  equivalent  first  order  system  in  an  appropriate  product  space,  we  can,  in  principle,  ap¬ 
ply  the  abstract  approximation  theory  developed  in  [8],  However,  since  the  resulting  first 
order  dynamics  are  not  in  general  described  by  a  strongly  monotone  operator,  the  results 
we  present  below  cannot  be  obtained  via  a  direct  application  of  our  first  order  theory. 
Only  the  essential  ideas  underlying  the  general  abstract  Banach  space  theory  involving 
m-accretive  operators  given  in  Section  2  of  [8]  are  directly  applicable. 


In  Section  2  below  we  review  the  abstract  existence,  uniqueness,  regularity,  and  approx¬ 
imation  results  for  nonlinear  evolution  equations  in  Banach  space  with  dynamics  governed 
by  m-accretive  operators.  In  the  third  section  we  consider  inverse  problems  for  abstract 
wave  equations  with  nonlinear  damping.  We  reformulate  these  second  order  equations  as 
equivalent  first  order  systems,  and  discuss  existence  and  regularity  of  solutions.  The  ap¬ 
proximation  theory  and  convergence  results  are  presented  in  Section  4.  In  the  concluding 
section  we  illustrate  the  application  of  our  theory  with  an  example  involving  estimation  of 
nonlinear  velocity  dependent  damping  in  a  one  dimensional  wave  equation. 


2.  Abstract  Nonlinear  Evolution  Equations  in  Banach  Spaces 


Let  Xo  be  a  Banach  space  with  norm  |-|0.  Let  T  >  0  and  suppose  that  Fo  G 
Li(0,T;Xo),  Let  A0  :  X o  — ►  2X°  be  an,  in  general,  nonlinear  set-valued  operator  with 
Dom(Ao)  =  {x  €  Xq  :  Aox  ^  0}.  We  assume  that  for  some  u  €  R  the  operator  Ao  +  u >1 
is  m-accretive.  That  is,  that  (i)  |ii  -  x2\Q  <  |(1  +  Aw)(ii  -  x2)  +  \{y\  —  j/2)|o  for  every 
x\,xi  €  Dom(>lo),yi  €  A0x  1,  y2  G  Aox2,  and  A  >  0,  and  (ii)  TZ(I  +  A (A0  +  u>I))  = 
U*eDom(/t0)  +  ^(Ao  +wl))x  =  Xo  for  some  A  >  0.  It  is  well  known  (see  [10]}  that  the 
operator  Ao  +wl  being  m-accretive  implies  that  for  each  A  >  0  the  resolvent  of  Ao  +w/  at 
A,  J(A;  Ao  4-u;/) :  X  — ►  X,  given  by  <7(A;  Ao  +wJ)  =  (I  +  A(^o  is  a  single  valued, 

everywhere  defined,  nonexpansive  (nonlinear)  operator  on  Xq. 
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We  consider  the  quosi-autonomous  initial  value  problem 

(2.1)  x(t)  +  A0x(t)  3  Fo(<),  0  <  t  <  T 

(2.2)  *(0)  = 

with  x°  G  Dom(A0).  By  a  (strong)  solution  to  the  initial  value  problem  (2.1),  (2.2)  we 
shall  mean  a  strongly  continuous  function  x  :  [0,  T]  — ♦  Xo  which  is  absolutely  continuous 
on  compact  subintervals  of  (0,T),  differentiable  almost  everywhere  on  (0 ,T),  and  satisfies 
F0(t)  -  x(t)  G  Aox(t)  for  almost  every  t  G  (0 ,T)  and  x(0)  =  a:0. 

A  two  parameter  family  of  nonlinear  operators,  {Uo(t}s)  :  0  <  s  <  t  <  T},  defined  on  a 
subset  ft  of  Xq  is  said  to  be  a  nonlinear  evolution  system  on  ft  if  it  satisfies  (i) 

Uo(t,  s)Uo(s,  r)<p  =  Uo(t,  r)p  for  every  ip  €  ft  and  r,  s,  t  with  0  <r<s<t<T,  and  (ii)  the 
mapping  (<,  s)  -»  C/0(f ,  ^)v?  is  continuous  from  the  triangle  A  =  {(r,  a)  :  0  <  a  <  r  <  T}  C 
R2  into  Xo  for  each  <p  €  ft. 

It  can  be  shown  (see  [8],  [12],  and  [14])  that  with  A0,F0,  and  x0  as  they  have  been 
defined  above,  there  exists  a  unique  nonlinear  evolution  system  {Uo(t,s)  :  0  <  s  <  t  <  T) 
on  Dom(Ao)  which  satisfies 

(i)  \Uo(t,  s)p  -  U0(t,  s)xl> |o  <  e | p  -  ^|0,  for  all  <p,  €  Dom(A0)  and  (t ,  s)  G  A, 

(ii)  \U0(s  + 1,  s)<p  -  U0(r  + 1,  r)v?|0  <2  /  \F0  {r  +  s)  -  F0  (r  +  r)|0  dr,  for  all 

o 

tp  G  Dom(Ao)  and  all  t  >  0  such  that  s  +  i,r  +  t  <T, 

(iii)  If  the  initial  value  problem  (2.1),  (2.2)  has  a  strong  solution  x ,  then  x(t)  =  t/0(M):r(s) 
for  all  (f,s)  G  A. 

The  strongly  continuous  function  x  defined  by  x(t)  =  Uo(t^O)x0,t  G  [0, T],  is  sometimes 
referred  to  as  the  unique  mild  or  generalized  solution  to  (2.1),  (2.2).  It  is  immediately 
clear  from  (iii)  above  that  when  the  initial  value  problem  (2.1),  (2.2)  has  a  strong  solution, 
it  coincides  with  the  unique  mild  solution  defined  in  terms  of  the  nonlinear  evolution 
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system  {£/<)(*,'*)  •  0  <  s  <  t  <  T}  on  Dom(/lo).  Theorem  2.1  below  is  the  fundamental 
abstract  approximation  result  upon  which  the  approximation  and  convergence  theory  for 
second  order  systems  to  be  developed  below  will  be  based.  It  is  similar  in  spirit  to  the 
approximation  results  for  nonlinear  evolution  systems  given  in  [13]  and  [15]. 

For  Banach  spaces  X  and  Y  we  denote  the  Banach  space  of  continuous  linear  operators 
from  X  into  Y  by  C(X,  Y). 

Theorem  2.1  For  each  n  G  Z+  =  {1,2, 3, . . . },  let  Xn  be  a  Banach  space  with  norm 
|  •  |„.  Let  7r„  G  C(Xo,Xn)  with  <  v  [<^>|0  where  u  is  a  positive  constant  independent 

of  n.  Let  Fn  G  Li(0,T;Xn)  and  let  An  :  Xn  -*  2Xn  be  an  operator  on  Dom(,4n)  C  Xn 
with  An  +ul  m-accretivc.  Suppose  that  there  exists  a  function  g  G  Li(0,T)  for  which 
|Fn(/)|n  <  </(t),  for  almost  every  t  G  (0 ,  T)  and  that 

(i)  for  each  tp  G  Dom(j4o)  there  exists  with  <pn  G  Dom(/ln)  such  that 

lim  | <pn  -  7rn^|n  =  0, 

ft—*  oo 

(ii)  lim  \Fn(t)  -  7rnF0(<)|n  =  0,  a.e.  t  G  (0,T), 

n— *oo 

(iii)  for  some  A  >  0,  we  have  lim  \J  (A;  An  +  ul)  <pn  —  nnJ  (A;  Aq  +  wl)  <p\n  —  0  whenever 

n  “■♦oo 

€  Xn  and  ip  G  Xo  with  lim  | pn  -  =  0-  Then  if  {Un  (t,s)  :  0  <  s  <  t  <  T}  is 

n— »oo 

the  evolution  system  on  Dom(A„)  generated  by  An  and  F„,  we  have 

(2.3)  lim  \Un(t,s)<Pn  -  irnU0(t,s)tp\  =0 

n— »oo 

uniformly  in  s  and  t  for  ( t,s )  G  A,  whenever  ip  G  Dom(>io)  and  <pn  €  Dom(>ln), 

with  lim  \<pn  -  7rny>|n  =  0. 
n  "  >  oo 

The  proof  of  Theorem  2.1,  which  we  omit,  is  essentially  the  same  as  the  proof  of  Theorem 
2.2  in  [8].  Some  modification,  albeit  relatively  simple  and  straightforward,  is  necessary 
however  due  to  the  fact  that  in  our  application  of  Theorem  2.1  below  we  cannot  choose  the 
X^b  as  subspaces  of  Xo.  In  the  proof  this  is  handled  by  simply  expressing  all  convergence 
conditions  and  results  in  terms  of  the  mappings  n„.  (Compare  the  statements  of  Theorem 
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2.1  above  and  Theorem  2.2  in  [8]).  We  note  that  the  convergence  result  (2.3)  is  closely 
related  to  the  notion  of  factor  convergence  os  given  in  [10],  We  remark  that  this  situation 
( Xn  X0 )  arises  frequently  in  parameter  estimation  problems  and  the  reader  can  find  a 
more  complete  discussion  with  examples  along  with  a  linear  semigroup  version  of  the  above 
approximation  theorem  in  [5,  Chapter  II].  Whether  or  not  (2.3)  implies  that  Un(t,s)<pn 
and  Uo(t,s)tp  are  close  in  any  reasonable  or  useful  sense  depends  of  course  upon  how  the 
spaces  Afn,  the  norms  |  •  jn,  and  the  mappings  7rn  are  chosen.  In  our  treatment  below 
we  shall  be  able  to  apply  the  following  corollary,  which  follows  from  Theorem  2.1  via  a 
straightforward  application  of  the  triangle  inequality. 

Corollary  2.1  Suppose  that  there  exists  a  Banach  space  X  with  norm  |  •  \x  for  which  Xo 
is  set- equivalent  to  X ,  for  each  n  =  1,2,...  ,  X„  is  a  subspace  of  X,  .and  all  of  the  norms 
I  .  |n,  n  =  0, 1, 2, . . .  are  uniformly  equivalent  to  |  •  |x.  Suppose  further  that  lim  TCny>  = 

n — >oo 

£  X.  Then  conditions  (i)  -  (iii)  in  the  statement  of  Theorem  2.1  can  be  replaced  by 

(i)  Dom(Ao)  C  ft.^‘00Dom(An)  ,  where  Li  denotes  the  lower  limit  of  the  sequence  of  sets, 
{Dom(An)}  (see  [11],  p.  335), 

(ii)  lim  Fn(t)  =  Fo(t )  in  X  for  a.e.  t  £  (0,  T),  and 

TI—+OO 

(iii)  for  some  A  >  0,  lim  J  (A;  An  +  wl)  tpn  =  J  (A;  An  +  u>I)  </?  in  X  for  each  <p  £  X  when- 

n  — »oo 

ever  <pn  £  Xn  with  lim  <pn  =  tp  in  X.  Moreover,  the  conclusion  (2.3)  can  be  replaced 

n— *oo 

by 

lim  Un(t,s)<pn  =  U0{t,s)ip 

n— *oo 

in  X ,  uniformly  in  t  and  s  for  (t,  s )  €  A,  for  each  ip  €  Dom(>lo)  and  <pn  £  Dom(  An),  n  = 
1,2,...  with  lim  pn  =  <p  in  X . 

fl— *00 

3.  The  Identification  of  Second  Order  Systems 

Although  the  identification  of  nonlinear  dissipation  mechanisms  in  infinite  dimensional 
second  order  (mechanical)  systems  is  of  primary  interest  to  us  here,  we  treat  the  somewhat 
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more  general  problem  which  includes  the  simultaneous  estimation  of  linear  stiffness.  It  will 
become  clear  from  our  discussions  below  that  an  approach  similar  to  the  one  we  use  here 
could  be  used  to  develop  techniques  for  the  estimation  of  nonlinear  stiffness  in  the  presence 
of  linear  damping.  The  extension  of  our  general  approach  to  handle  inverse  problems 
involving  the  simultaneous  identification  of  both  nonlinear  stiffness  and  damping  is  at 
present  still  under  investigation.  Wc  hope  to  be  able  to  report  on  results  in  this  direction 
in  the  not  too  distant  future. 

Let  H  be  a  real  Hilbert  space  with  inner  product  <  v  >  and  corresponding  norm  |  •  |. 
Let  V  be  a  reflexive  Banach  space  with  norm  denoted  by  ||  •  ||,  and  assume  that  V  is 
densely  and  continuously  embedded  in  H .  If  we  let  V*  denote  the  dual  of  V  it  follows  that 
\r  <_+  H  <_♦  V*  with  the  embedding  of  H  in  V*  dense  and  continuous.  The  continuity  of 
the  embeddings  implies  the  existence  of  a  constant  /z  >  0  for  which  \ip\  <  fi  ||y>||  ,<p>  G  V, 
and  |M|.  <  where  ||-||,  denotes  the  uniform  operator  norm  on  V*.  We 

shall  also  use  <  •,  •  >  to  denote  the  usual  extension  of  the  H  inner  product  to  the  duality 
pairing  between  V  and  V*.  Let  Q  and  Z  be  metric  spaces  and  let  Q  be  a  fixed,  nonempty, 
compact  subset  of  Q. 

For  each  q  G  Q  let  the  operator  A(q)  €  C(V ,  V *)  satisfy  the  following  conditions: 

(Al)  (Symmetry)  For  all  €  F,  <  A(q)<p,ip  >=<  ^p,A(q)rp  >; 

(A2)  (Continuity)  For  each  9?  S  V  the  map  q  ->  A(q)(p  is  continuous  from  Q  C  Q  into  V*\ 

(A3)  (Equi-V-Coercivity)  There  exist  constants  ax  €  R  and  a0  >  0,  both  independent  of 
q  €  Q  for  which  <  A{q)tp, if  >  +<*1  M2  ^  IM|2»  ^or  V  e  ^  anc*  ?  e  Q' 

(A4)  (Equi-Boundedness)  The  operators  A(q)  ore  uniformly  bounded  in  q  for  q  G  Q]  that 
is,  there  exists  a  constant  a  >  0,  independent  of  q  G  Q,  for  which  ||<4(g)v’||*  <  a  |MI> 
for  all  tp  G  V. 
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In  addition,  for  each  q  G  Q  let  the  operator  B(q)  :  T)om(B(q))  C  V  — ♦  2'/*  satisfy  the 
following  conditions: 

(Bl)  (Domain)  Dom(5(tf))  =  Dom(S)  is  independent  of  q  for  q  G  Q,  and  0  €  Dom(B); 

(B2)  (Continuity)  For  each  y?  G  Dom(5),  the  map  q  — ♦  (B(q)tp  is  lower  semi-continuous  from 
Q  C  Q  into  2V'  (see  [17],  Section  18.i,  page  173); 

(B3)  (Maximal  Monotonicity)  For  all  (y>i,  tpi) ,  (^2^2)  6  Bq  =  {(y?,  tp)  G  V  x  V*  : 

y?  G  Dom  ( B),i/>  G  &(q)<p}  we  have  <  Vh  -  ^2, ¥>i  “  ^2  >  >0  with  not  properly 
contained  in  any  other  subset  of  V  x  V*  for  which  this  monotonicity  condition  holds; 

(B4)  (Equi-Boundedness)  The  operators  B(q)  map  V-bounded  subsets  of  Dom(B)  into  subsets 
of  V*  which  are  uniformly  V*-bounded  in  q  for  q  G  Q\  that  is,  if  5  is  a  V-bounded  subset 
of  Dom(5),  then  the  set  B(q)S  is  V*-bounded,  uniformly  in  q  for  q  G  Q. 

For  each  q  G  Q  let  u°(q)  G  V,ul  (q)  G  H  and  /(•;?)  6  Lx  (0 ,T;H)  where  T  >  0  is 
given  and  fixed.  We  shall  require  that  the  mappings  q  — »  u°(q)  £  V,q  u*( q )  G  H ,  and 
q  -*  /(<;?),  for  almost  every  t  G  (0 ,T),  are  continuous  from  Q  C  Q  into  V,H,  and  H 
respectively.  For  each  z  €  Z,  let  <j>  (•  ;  z)  be  a  continuous  mapping  from  C  (0,  T;  V  x  H) 
into  R+  and  consider  the  identification  problem  given  by: 

(ID)  Given  observations  z  G  Z,  determine  parameters  q  G  Q  which  minimize  the  performance 
index 

*(«)  =  *  ((«(*),  *(*));*) 

where  for  each  q  G  t?,u(g)  —  u(*  ;</)  is  the  mild  solution  to  the  second  order  initial 
value  problem 

(3.1)  «(<)  +  B(q)u(t)  +  A(q)u(t)  9  f(t;q) ,  0<t<T 

(3.2)  u(0)  =  «°(g),  u(0)  =  «■(«). 


8 


Following  Barbu  [10],  we  rewrite  the  second  order  initial  value  problem  (3.1),  (3.2)  as 
an  equivalent  first  order  system  in  a  product  space  and  then  apply  the  abstract  theory 
outlined  in  Section  2  above  to  make  the  notion  of  a  mild  solution  precise.  Let  X  •—  V  x  H 
be  the  Banach  space  endowed  with  the  norm  |  •  j  given  by 

Kv.WIx  =  (lMIJ  +  Ms)5  ■ 

In  addition,  for  each  q  €  Q  let  H(q)  be  the  Hilbert  space  which  is  set  equivalent  to  X  and 
has  as  inner  product  <  •,  •  >q  given  by 

(3.3)  <  (y?1,V>i),  (92,^2)  >,  =  <  A(q)<f 1,^2  >  +<*1  <  <pu  VI  >  +  <  01,t/»2  >  . 

Denote  the  corresponding  induced  norm  by  |  •  |?.  Conditions  (Al)  -  (A4)  above  guarantee 
that  (3.3)  indeed  defines  an  inner  product  on  X  and  that  the  Banach  spaces  {’ H(q ),  |- } 
are  uniformly  (in  q ,  for  q  G  Q)  norm  equivalent  to  the  space  {X,  |  •  |x}- 

For  each  q  G  Q  define  the  operator  A(q)  :Dom(A  (q))  C  'H(q)  -+  2*^  by 
’  A  (q)  (y?,0)  =  {-'!>,{A(q)<p  +  B(q)*l>}  D  H) , 
with  Dom(A(g))  =  {(y?,0)  €  V  x  V  :  0  G  Dom(0) ,  {A(q)<p  +  B  ( </)0 }  fl  H  7^  0}. 

Theorem  3.1  There  exists  an  u  G  R,  independent  of  q  G  Q,  for  which  the  operator 
A(q)  +  ujI  is  m-accretive. 

Proof.  Our  proof  is  analogous  to  the  one  given  by  Barbu  [10]  in  the  case  where  the 
space  and  operator  do  not  depend  upon  a  parameter.  We  first  show  that  for  u  G  R 
sufficiently  large  the  operator  A(q)  +u>I  is  monotone.  For  (<^»,  0) ,  (77,  G  Dom(A(^)) 
and  u  G  R  let  (0,0)  G  (A(q)  +  and  (77, 0)  G  (A  (q)  +  w/)  (77,  6).  It  follows  that 

A  A 

0  =  -0  +  u<p,  fj  —  -8  +  u tj,  0  =  A  (q)  y?  4-  0  4-  o>0  and  8  -  A  (q)  77  +  8  +  u>8  for  some 

AT  A» 

0  G  B  (q)  0  and  some  6  G  B(q)8.  Then,  using  the  monotonicity  of  the  operator  B(q) 
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(condition  (B3))  we  obtain 

<  (<£,V0  “  (*b  #)»  (?.  VO  “  (*?.  *)  >? 

=<  (-V’  +  0,*4(?)(9  -*?)  +  V*  -  0)>(<^  ~  f?>V»  -  d)  >q  +u\(<?-ri,ip  -0)\ ] 

-  -ai  <  ip  ~  0,tp  -  rj  >  +  <  tj>  -  0,rp  -  0  >  +U)\(<p  -  r],rp  -  6) |J 

>  -ct\  <  rp  -  0,  <f  -  r)  >  +u>  <  A  (q)  (v?  -  t]) ,  <p  -  rj  >  +ua\  |</?  -  t]\2  +  u  \xj>  -  0\ 2 

>  0 


if 


a;  >  uo  =  sup 


ai  <  u,u  >  ,  N  „ 

- s - r;  (u,  i>)  G 

<  A(q)u,u  >  +ai  |u|2  +  |i>|2 


We  note  that  u>o  exists  since 


X,(u,v)  ^  0,q  G 


ati  <  u,v  >  aifi  ||u||  |t>| 

<  A(q)u,v  >  +<*i  |u|2  +  M2  ~  a0  ||u||2  +  |u|2 

^  ?oi/i{H|a  +  M2}  ^  QlfX 

<*oN|2  +  M2  “  2min(a0>  1) 

If  u0  <  0,  we  set  u  =  0  and  henceforth  assume  without  loss  of  generality  that  u  >  0.  Since 
A(q )  +  is  monotone,  it  is  also  accretive.  Indeed,  letting  A  >  0  ,  we  see  that 

l(v>»  VO  “  (*? , 0)\*  <  |(v»,  VO  -  fa  8)\] 

+  A  <  (£,  fy  -  ( rj ,  0),  (y>,  V»)  -  (??,  0)  >q 

=<  (v>,  VO  -  (*7,  +  H(<P,  VO  -  (v,  #)},  (v>,  VO  ~  (*?,  0)  > q 
<  l(¥>>  V’)  -  (*7,  *)  +  M(<£,  ^ )  “  (»?,  0)}|?  |(p,  V»)  -  (17,  0)1,  ■ 

Thus  we  obtain 


l(v>,V0  -(*7,0)1,  <  l(^VO-(^0)  +  M(£,VO“(»7>0)}|?- 

To  demonstrate  m-accretivity,  we  let  A  >  0,  and  let  (77,  9)  G  be  fixed  but  arbitrary, 
and  show  that  (I  +  \(A(q)  +  u>I))(<p,il>)  3  (77,0)  for  some  (y?,V0  G  Dom(A(g)).  This 
statement  can  be  equivalently  written  as 

(3.4)  <p  +  +u)y>)  =  rj 
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(3.5) 


ip  4-  A  (A  (?)  tp  4-  B(q)rp  +  wi p)  3  6. 


Solving  for  q>  in  (3.4)  and  the  substituting  into  (3.5),  we  obtain 

(3.6)  =  (1  +  Aw)  1  (A  0  +  r?) 

(3.7)  (1  +  Aw)  rP  +  A2  (1  +  Aw)"1  A  (?)  0  +  XB  (?)  xP  3  6  -  A  (1  +  Aw)"1  A  (?)  77. 

Define  the  operator  T\  (?)  £  £(V,  V*)  by 

(3.8)  Ta  (?)  =  (1  +  Aw)  /  +  A2  (1  +  Aw)"1  A  (?) . 

It  follows  from  condition  (A3)  that  for  u  £  V 

(3.9)  <  T\  (?)  u,u  >  >  Tj  (A)  |u|2  +  r0(A)  ||u||2  , 

where  rj  (A)  —  (1  -f  Aw)  —  A2  (1  +  Aw)-1  «i  and  To  (A)  =  A2  (1  +  Aw)-1  oo-  F°r  A  >  0 
appropriately  chosen  (small  if  <*1  >  0;  any  A  >  0  if  c*i  <  0),  it  is  seen  that  ri(A)  >  0  and 
hence  that  T\  (?)  is  monotone  and  coercive  in  the  sense  of  Barbu  [10,  p.  34].  That  is, 

<  T\(q)u,u  >  /  ||u||  -+  00,  as  ||u||  00. 

Therefore,  applying  Barbu’s  [10]  Corollary  II. 1.3  with  his  X ,  B  and  A  chosen  as  X  = 
V,B  =  T\(q)  and  A  =  A £(?),  we  find  that  ll(T\(q)  +  A B(q))  =  V*.  Consequently  there 
exists  a  rp  £  Dom(5  (?))  for  which  (3.7)  holds.  If  we  then  obtain  if  from  (3.6),  it  follows 
that  71(1  +  A(A(?)  +w/))  =  H  (?)  and  the  theorem  is  proved. 

For  each  ?  £  Q  define  F  €  L\  (0 ,T;7f  (?))  by  F(i\q)  =  (0, /(<;?)),  a.e.  i  £  (0,T),  and 
set  x°  (?)  =  (u°(?),u1  (?))  £  H(q).  Theorem  3.1  and  the  discussions  of  Section  2  yield 
that  A(q)  and  F(*;  ?)  generate  a  nonlinear  evolution  system  {17  (f,  s;  ?) :  0  <  s  <  t  <  T}  on 
Dom(A(?)).  Henceforth  we  shall  assume  that  x°(?)  £  Dom(A(?))  for  each  ?  £  Q,  and  by  a 
mild  solution,  it(?),  to  the  second  order  initial  value  problem  (3.1),  (3.2)  we  shall  mean  the 
V-continuous  (recall  the  uniform  norm  equivalence  of  X  and  H(q))  function  u(*;  ?)  given  by 
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the  first  component  of  the  H(q)  (or  X)  continuous  function  x  (<;  q)  =  U  (t,  0;  q)  x°  (q) ,  t  € 
[0,T].  We  shall  take  u('\q)  to  be  the  /f-continuous  second  component  of  x(-;q). 


At  this  point  some  remarks  regarding  the  closure  of  the  set  Dom(A(2))  are  in  order.  We 
show  that  under  a  slightly  more  restrictive  boundedness  condition  on  the  operators  B(q) 
than  condition  (B4),  the  operator  A(q)  is  in  fact  densely  defined;  that  is,  that  Dom(A(<7))  = 
H{q). 

Theorem  3.2  If  condition  (-B4)  is  replaced  by  the  stronger  condition  that  the  operators 
B(q)  map  if-bounded  subsets  of  Dom(B)  into  subsets  of  V *  which  are  uniformly  V*- 
bounded  in  q  for  q  6  Q,  then  Dom(A(g))  =  T~i(q)- 

Proof:  Let  (y>,i/>)  £  'H(q)  and  set 

ten,  V>n)  =  J  ihA(<l)+  uI)  te>  VO  €  Dom  (A  (q))  , 
for  n  =  1, 2, - As  in  the  proof  of  Theorem  3.1,  we  find 

(3.10)  (pn  =  (l  +  n-1u;)  1  (n"Vn  +  </>) 

(3.11)  (Tn-i  (?)  +  n_1S  (?))  3  !/>  -  n"1  (l  +  n"'w)  1  A{q)p> 

where  the  operator  Tn- 1  (q)  e  £(V,  V*)  is  given  by  (3.8)  with  A  =  n-1.  Recalling  (3.9), 
we  have  for  u  €  V 


<  Tn-\  (q)u,u  >  >  T\  (n  ^luf+Tote  a)  |jw||2  , 

where  T\  ( n~l )  =  (l  +  n-1w)  —  n-2  (l  +  n-1u ;)  1  Oi  >  0  for  n  sufficiently  large,  and 

=  n”2(l  +  n_1u;)  1  aQ.  Note  that  lim  rx  (n~l)  =  1  and  lim  r0  (n_1)  =  0. 

v  '  v  n— *oo  x  '  n— *oo  x  ' 

Now  from  condition  (B3)  we  find  that  for  some  i j>n  £  B  (q)  ipn  and  any  6  €  B  (q)  (0) 
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(recall  that  0  6  Dom(8)  by  condition  (Bl)), 

r\  (n_1)  |^„|2  +r0  (n"1)  ||i/>„||2  <  <  Tn-i  (q)ipn,rpn  > 

<  <  Tn-i  {q)ipn^n  >  +n_1  <  rpn  -  0,ipn  > 

=<  ip  —  n-1  (l  +  n-1w)  1  A(q)  <p  —  n~lQ,  ipn  > 

<  M HM  +  n"1  (1  +  n"1^)"1  |j A  (?)  y\\.  ||0„||  +  n"1  ||0||*  (|^«| 

<  M2  +  ^  M2  + 


+ 


^ll^||2  +  i^||«||5  +  ^P 


ll<M2. 


Therefore 


iw2 + iiwi2  <  w2  +  i  im(«mi  2  +  ^m2.- 


It  follows  that  {ipn}  is  uniformly  bounded  in  H  and  that  {n  1  \\ipn\\ }  is  also  uniformly 
bounded.  Prom  (3.11)  we  obtain 


ipn-ip  G  -n~luipn~n~2  (l  +  n  1oo)  1  A(q)ipn~n  1  (l  +  n-1^)  1  A(q)ip-n  lB(q)ipn- 


Therefore  for  any  ipn  €  B  (q)  ipn,  we  have 

II V’n  -V’lL  <  n~lu  HV'nll*  +  n~2  (l  +  n_1u;)  1  ||.4(g)  V’nll, 

+  n_1  (1  +  n-1u/)  1  |M(g)v?||*  +  n'1  \\xpn\\* 

<  n~'un\ipn\  +  n_2a  HV’nll  +n_1a||y>||  +  n_1 1|^„||*  . 

This  estimate  together  with  the  uniform  boundedness  of  IV’nl  and  n~l  \\ipn ||  and  the  hypoth¬ 
esis  of  the  theorem  imply  that  lim  ipn  =  ip  in  V*.  Thus  using  the  uniform  boundedness 

n—*oo 

of  {IV’nl}  and  the  density  of  V  in  H,  we  obtain  ipn  — >  xp  weakly  in  H  as  n  — ►  oo. 


From  (3.10)  we  find  that 


M  <n  1  |tf„|  +  M 


and 

ll^nll  +  IMI 
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and  therefore  that  the  {y?,,}  are  uniformly  bounded  in  both  H  and  V ,  Since  (l  +  n_,u>)  ipn 
=  n~ltpn  +  <p,  we  have 

| <fin  ~  <p\  <  ”_1  iV'nl  +  n~l W  \tpn\ 

<  (n_1  +  n"2u>)  |t/>n|  +  n_1u;  \p\  — ♦  0 

as  n  — »  oo.  Thus  lim  y?n  =  y?  in  if.  Using  the  density  of  H  in  V*  and  the  uniform 

n—*oo 

V-boundedness  of  the  ipn,  it  is  not  difficult  to  argue  that  < pn  — ♦  tp,  weakly  in  V.  Therefore 
(y?n,  V>„)  — ►  (y?,t/>)  weakly  in  X  and  thus  the  weak  closure  of  Dom  (A(</))  is  all  of  X.  Since 
A(q )  4-wI  is  m-accretive  by  Theorem  3.1,  the  strong  closure  of  Dom  ( A(q ))  is  convex.  It 
follows,  therefore,  that  Dom(A(g))  =  'H(q)  and  the  theorem  is  proved. 

When  we  consider  examples  in  Section  5  below,  we  shall  show  that  there  is  an  im¬ 
portant  class  of  nonlinear  models  for  dissipation  which  satisfy  the  boundedness  condition 
hypothesized  in  the  statement  of  Theorem  3.2. 

Finally,  we  note  that  another  class  of  operators  B(q)  which  lead  to  the  operator  A(q) 
being  densely  defined  is  the  one  which  includes  operators  B(q)  G  C(V,  V*)  for  which  the 
set  {p  €  V  :  B(q)p  G  H }  is  dense  in  H.  For  then  condition  (A3)  implies  that  {p  G  V : 
A(q)<p  G  H]  is  dense  in  V  (see  [22],  Theorem  III.2.B)  and  consequently  Dom(A(5))  = 
{(y?,V0  G  V  xV:  A(q)y>  +  B(q)ip  G  H}  is  dense  in  7 i(q).  In  particular,  this  will  be  the 
case  for  the  well  known  linear  Kelvin- Voigt  (i.e.,  B{q)  ~  -4(9))  and  the  so-called  structural 
(i.e.  B(q)  ~  A(q)?)  viscoelastic  damping  models. 

4.  Approximation  Theory  and  Convergence  Results 

We  apply  Ritz-Galerkin  finite  element  techniques  to  discretize  the  state  equation.  For 
each  n  =  1,2,...  let  Hn  be  a  finite  dimensional  subspace  of  H  with  Hn  C  V.  Let 
Pn  :  H  Hn  denote  the  orthogonal  projection  of  H  onto  Hn  computed  with  respect  to 
the  H  inner  product.  We  assume  that  PnDom(5)  C  Dom(B)  and  that  the  approximation 
condition 

(P)  lim  \\Pn<p  -  y>||  =  0,  ¥>  €  V, 

n— *oo 
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is  satisfied.  We  note  that  condition  (P)  implies  that  lim  | Pn<p  —  <p\  —  0,  v3  €  H  as  well. 

fj— *00 

For  each  n  =  1,2,...  and  q  £  Q)  let  An(q)  €  C(Hn)  =  C(Hn,Hn)  and  Bn{q)  :Dom 
(Bn)  C  Hn  —*  2Hn  denote  the  usual  Galerkin  approximations  to  the  operators  A(q)  and 
B(q)  respectively.  More  precisely,  for  </?„  £  Hn  let  An{q)<Pn  =  V’n  where  ip„  is  the  unique 
element  in  Hn  guaranteed  to  exist  by  the  Riesz  representation  theorem  satisfying 
<  A(q)<pn,Xn  >  =  <  V>n,Xn  >,X«  £  H„.  Similarly,  for  <pn  £Dom(Bn)  =  Dom(5)n/f„  f 
0,  let  Bn{q)q>n  =  {xpn  £  Hn  :<  ij>, Xn  >=<  'I’n'Xn  >,Xn  €  Hn,  for  some  ip  £  B(q)ipn}  ■ 
We  set  /„(•;<?)  =  Pnf{-\q)  €  L,  (0 ,T;  ffn),  u°„(9)  =  P„u°(9)  and  u»(g)  =  P^1^). 

Our  approximation  framework  is  based  upon  the  following  sequence  of  approximating 
identification  problems. 

(IDn)  For  the  observations  z  £  Z  given  in  problem  ( ID ),  determine  parameters  qn  €  Q  which 

minimize  the  functional 

#n(<?)  =  <Kun(q)yun(q)\z) 

where  for  each  q  £  Q  un(q)  =  u„(-;g)  is  the  mild  solution  to  the  second  order  initial 

value  problem  in  Hn  given  by 

(4.1)  un(t)  +  Bn(q)un(t)  +  An{q)un(t)3fn(t-,q),  0  <t<T, 

(4.2)  u„(0)  =  u°„(j),  u„(0)  =  «  i(«). 


Once  again  to  establish  the  existence  of  a  unique  mild  solution  for  each  q  £  Q,  and  to 
develop  a  convergence  theory,  we  reformulate  the  initial  value  problem  (4.1),  (4.2)  as  an 
equivalent  first  order  vector  system  in  an  appropriate  product  space,  and  then  apply  the 
abstract  existence  and  approximation  results  outlined  in  Section  2. 

For  each  n  =  1,2, . . .  let  Xn  =  Hn  x  Hn  and  for  each  q  £  Q  let  Hn{q)  denote  the  space 
Xn  considered  as  a  subspace  of  the  Hilbert  space  H(q).  Define  the  operator  An{q )  :  Dorn 
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(An(q))  c  Hn(q)  -  by 


An(q)  (^n,  V>n)  “  (-^n,X(?K  +  Bn(q)^n) 

for  each  (y>n,0n)  €  Dom(.An(tf))  =  if„x  Dom (tfn).  Set  F„  (t;?)  =  (0,/n  (t;g))  for  almost 
every  t  €  (0,T)  and  let  x°n(q)  =  (9), «),(?)).  Note  that  Fn(-‘,q )  €  L\  (0 ,T-,Hn(q))  and 
that  the  assumptions  that  x°(q)  6  Dom(i4(g))  and  Pn  Dom(5)  C  Dom  ( B )  imply  that 
x°n(q)  €  Dom(A(q)). 

With  the  above  definitions  of  An(q)  and  Bn(q ),  it  is  not  difficult  to  argue  that  An(q)+u>I 
is  m- accretive  in  'Hn(q).  (The  proof  is  analogous  to  the  proof  of  Theorem  3.1.)  It  follows 
that  for  each  n  =  1, 2, . . .  and  each  q  €  Q,  An(q )  and  F„(-;  q)  generate  a  nonlinear  evolution 
system,  {Un  (t,  s;  q)  :  0  <  s  <  t  <  T)  on  Dom(.4(g)).  By  a  mild  solution  un(q )  =  ««(•;?) 
to  the  initial  value  problem  (4.1),  (4.2)  we  shall  mean  the  V-continuous  first  component  of 
in  (-;q)  =  U„  (-,0 -,q)x^(q).  The  if-continuous  second  component  of  x„(-;<?)  will  be  taken 
to  be  u„  (•;$). 

The  primary  result  of  this  paper  is  given  in  the  following  theorem. 

Theorem  4.1.  For  each  n  =  1,2,...  the  problem  ( IDn )  admits  a  solution  qn  €  Q. 
Moreover,  the  sequence  {gn}^Li  has  a  convergent  subsequence  {qnk  }^Li  with  lim  qnk  = 

k — >oo 

q  €  Q,  where  q  is  a  solution  to  problem  (ID). 

Proof:  In  light  of  the  compactness  of  Q  and  the  continuity  of  <f>,  the  existence  of  a  solution 

qn  €  Q  to  problem  ( IDn )  can  be  established  by  demonstrating  that  the  mapping  q  — * 
xn(-  ;  q )  is  continuous  from  Q  C  Q  into  C(0,  T;  X).  Furthermore,  as  we  have  demonstrated 
in  several  of  our  earlier  papers  (see,  for  example,  [6],  [8]),  showing  that  q  =  lim  qnic  is 

k—>oo 

indeed  a  solution  to  problem  (ID)  requires  arguing  that  for  any  convergent  sequence 
{tfn}^i  with  qn  €  Q  and  lim  qn  =  q0  e  Q  we  have 

n — *oo 

(4.3)  Jiimo|xn(t;9n)-x(t;g0)|x  =0 
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uniformly  in  t  for  t  G  [0,T].  Both  the  continuous  dependence  and  the  convergence  (4.3) 
can  be  demonstrated  via  an  application  of  Corollary  2.1.  We  Bhall  establish  (4.3)  here; 
continuous  dependence  can  be  argued  in  an  analogous  manner. 

Let  {tfnirLi  C  Q  with  lim  qn  —  qo  in  C*  The  compactness  of  Q  of  course  implies  that 

fl  — *oo 

<7o  6  Q.  In  applying  Theorem  2.1,  we  take  X0  =  'H(qo),  Xn  =  ft„(<?n),  Mr,  =  H?n  *n  = 
0,1,2,...  ,  and  define  7r„  G  £ (W(<7o) ,  ft* (?«))  by  7r„(<^,0)  =  (P„v?, Pnt/») , (y>, r/»)  G 
H(qo).  Note  that  condition  (P)  implies  that  the  7rn  are  uniformly  bounded  and  converge 
strongly  in  X  to  the  identity.  We  set  Fn  =  Fn  (■;qn)  ,F0  =  F(-;q0),A„  =  An(qn),A  0  - 
A(q0),  and  will  use  Corollary  2.1  to  establish  that 

lim  |Un(i,0;gn)a:°  (qn)  -  U  (i,0\qo)x°  (q0)\x  =  0 

uniformly  in  rf,  for  t  G  [0,  T].  Under  the  general  assumptions  that  we  have  made  above,  all 
of  the  hypotheses  of  Corollary  2.1  can  be  verified  immediately  with  the  exception  of  the 
resolvent  convergence  condition  (ii).  We  establish  it  here. 

Let  A  >  0  be  fixed,  let  (</5,i/>)  €  H(qo),  (Vmi’n)  €  %„(<?„),  and  assume  that 
lim  Kv^mV’n)  —  (<P,*I>) \x  —  0.  It  then  follows  that  lim  </?„  =  y?  in  V  and  lim  ipn  —  ip  in 

n— *oo  n— *oo  n — >oo 

H.  Setting  ( r/,9 )  =  J(\]  A(q0)  +  ul)(<f>,ij>)  G  Dom  (A(g0))  and 
{Vn,0n)  =  j  (A;  An  ( qn )  +  wl)  (<^>n,  4'n)  G  Dom(A„  (qn))  we  obtain 

(4.4)  77  =  (1  +  Aw)-1  (X6  +  ip) 

(4.5)  ( Tx  (qo)  +  \B  (q0))  9  3  xj>  -  A  (1  +  Aw)"1  A(q0)<p 

(4.6)  qn  =  (1  +  Aw)-1  (\9n  +  </>„) 

(4.7)  (Pn.A  (<?n)  +  A5n  (qn))  @n  ^  iftn  ~~  A  (1  +  Aw)  An(qn)ipn, 
where  T\(q)  is  given  by  (3.8),  and 

Tn, A  (q)  =  (1  +  Aw)  /  +  A2  (1  +  Aw)"1  An  (q)  G  C(Hn). 
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Recalling  (3.9),  we  choose  A  >  0  bo  that  r\(\)  >  0.  Wc  first  show  that  the  9n  6  H„ 
are  uniformly  V-bounded  in  n.  From  (3.9),  conditions  (A4),  (B3),  and  (B4),  and  the 
definitions  of  the  operators  A»(q)  and  B„(t j)  we  find  for  some  §n  6  Bn(qn)9n  and  any 

(n6B(?n)(0) 

ri(A)|0„|2  +  ro(A)||0„||2  <  <  Tx > 


—  TntX  (7n)  0ni  0fi  A 

^  <  Pnt\  (?n)  9 n ,  0 n  >  4" A  <  9n  —  £n>9n  > 


—  <  “  A  (1  +  Xu>)  A  ( qn )  “  A£n,  0n  > 

<  IV’nl  |0„|  +  A(1  +  Aa))”1  a  ||pR||  ||*n||  +  A/?o  ||0„|| 

<  _i_  i^.i2  +  \ej  +  £!il±MJ.||9n|| 


+  ^  ll«nf  +  (1  +  ^)g,2  +  ||0„||2 

4  Of0  4 


where  /?o  =  5up  1 ||£||*  :  £  €  i?  (9)  (0) ,  q  6  Q  j.  It  follows  that 


n(A) 


l^n|2  + 


r0  (A) 


m,2  <  ^  i<m2 + M2  + 


2  ■  2  -  2n  (A)  '  '  o0  .  o„ 

The  fact  that  lim  ||<pn  —  <p||  =  0  and  lim  \^n  —  ^1=0  implies  that  ||<pn||  and  \ipn\  are 

n~*oo  n—*oo 

uniformly  bounded,  and  consequently  that  the  9n  are  uniformly  V-bounded. 


Once  again  from  (3.9),  for  any  9n  (qn)  €  B  (qn)  9 n,  9  (qn)  €  B(qn)9,  and0(go)  €  B(q0)  9 
we  have 

T\  (A)  \e„  -  e\ 2  +  r0  (A)  ||#„  -  0||2  <  <  Tx  (?„)  {0„  -  0) ,0„  -  0  > 

<  <Tx  (?„)  (0„  -  0} ,  0„  -  0  >  +A  <  0„  ({„)  -  0  ( gn ) ,  0„  -  0  > 

=  ^  Trt(A  (?n)  9n  —  T\  (qo)  9 ,  0n  —  -Pn0  >  4*  A  <  0„  (9n)  ~  0 (90) )  9n  —  Pn9  > 

4-  <  T\  (qn)  9n  —  T\  (qo)  9 ,  Pn9  —  9  >  +A  <  0n  (<?„)  —  0  (90) ,  Pn9  —  9  > 

+  <  TA(9o)0 -TA(9„)0,0n  -0  >  +A  <  0(?o)  -  0(?n),0n  -  0  >  . 

Now  from  (4.7)  there  exists  a  9„  €  Bn(qn)9n  for  which 

Tn, A  (9n)0n  4"  A0„  =  l/>n  -  A  (1  +  Au>)_1  ^n^rO^n- 
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Choose  dn(g„)  to  he  any  element  in  B(q„)9n  for  which 


~  <  *  X  n  Xn€  Hn • 

That  such  a  9n(q„)  €  B(qn)8n  exists  is  “guaranteed"  by  the  definition  of  Bn{qn)-  Also, 
recalling  (4.5),  we  choose  9 (go)  £  B  (go)  8  so  that  T\  (go)  8  4  X9  (go)  = 
rp  —  A(1  4  Aw)”1  -4 (go) <£•  Then 

T\  (A)  \9n  —  9\2  4  To  (A)  j|0n  —  8\\2 

<  <  T„,a  (g»»)  8n  4  A 9n  —  {T\  (go)  8  4-  A 8 (go ) } » 9n  —  Pn9  > 

4  <  Tx  (gn)  8n  -  Tx  (g0)  9 ,  Pn8  -  8  >  +  A  <  9n  (gn)  -  8  (g0) ,  P„0  -  0  > 

4  <  Tx  (go)  0  —  Ta  (g„)Mn  “  *  >  +A  <  0(go)  -  0(g„)  A  -  9  > 

—  <  Ipn  ~  0,  $n  —  Pn9  >  4  A  (1  4  Aw)  1  <  .4  (go)  V’  —  «An  (?n)  <Pn)9n  —  Pn9  > 

+  <TX  (gn)  -  Ta  (go)  0,  Pn<?  -  9  >  4  A  <  8n  (gn)  -  9  (go)  ,Pn9-8> 

4-  <  Ta  (go) 9  —  Ta (gn) 9,9n  —  8  >  4  A  <  8 (go)  —  6* (gn) i8n  —  9  > 

=  <  0n  --  0, 0„  -  Pn9  >  +A  (1  -f  Aw)"1  <  .4  (g0)  ^  -  -4  (g„)  y?,  0n  -  Pn9  > 

4  A  (1  4  Aw)  <  A  (g«)  (y?n  —  ^p)  ?  9n  —  Pn9  >  4  (1  4  Aw)  <  9n  —  9,  Pn8  —  8  > 

4  A2  (1  4  Aw)"1  <  A  (gn)  $n  -  *4  (g0)  0,  Pn0  -  0  >  4A  <  0„  (gn)  -  9  (g0)  ,Pn9-8> 

4  A2  (1  4  Aw)"1  <  (A (g0) 0  -  A (g„) 0, 0„  -  9  >  4A  <  9 (g0)  -  0 (gn) A  -  8  > 


<  | i>n  -  011 9n  -  Pn8 1  4  A  (1  4  Aw)  1  ||4 (g0)  y?  -  A (gn)  y>||.  ||0n  -  Pn*|| 

4  A  (1  4  Aw)"1  M(gn)  (?«  -  <f)L  I A  ~  Pn9\\  4  (1  4  Aw)  \8n  -  8\  \Pn9  -  9\ 

4  A2  (1  4  Aw)"1  \\A(qn)8n  -  A(q0)9\i  \\Pn8  -  9\\  4  A|| 0„  (gn)  -  8 (q0)\\.\\Pn8  -  8\\ 
4  A2  (1  4  Aw)"1  \\A(qo)9-A(qn)9\l  |A  -  9\\  4  A||0(go)  -  9{qn)  ||.  || 8n  -  9\\ 
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<  {I*.  I  +  1*1)  W.  -<P\  +  A(1  +  Aw)"1  {||#„||  +  K  ||«|||  \\Mio)v  -  A(<l„)v\\. 

+  A(1  +  Au.)'1  or  {||#„||  +  1<  ||#||)  ||y>B  -  v»||  +  (1  +  Au>)  {|#„|  +  |#|)  \PJ  -  #| 

+  AJ  (1  +  Au.)-'  O  {||«„||  +  ||#||)  ||P„#  -  #||  +  A{||#„  (#„)||.  +  ||#(«o)||.}  || P.t  ~  #11 
+  A5  (1  +  Aar)'1  {||#„||  +  ||#||}  |M(!o)«-  a((,„)#||.  +  A  {||#„||  +  ||#||)  ||#(,o)  -  #(,„)||. 

where  K  is  the  uniform  bound  on  the  operators  Pn  6  C{V,V)  guaranteed  to  exist  by 
condition  (P).  By  condition  (B2)  we  can  choose  6(qn)  €  B{<ln)Q  so  that  lim  9(qn)  =  $(90) 

n— ♦  00 

in  V*.  Thus,  the  fact  that  lim  V>n  =  »n  H,  lim  tpn  =  V  in  V ,  the  uniform  boundedness 

n— *  00  n— *oo 

of  ||0„||  (and  therefore  |0„]  as  well),  the  continuity  condition  (A2),  and  the  boundedness 
condition  (B4)  imply  that  ^i^^n  =  9  in  both  H  and  V\  From  (4.4)  and  (4.6)  it  then 
immediately  follows  that  lim  n„  =  n  in  V  and  the  theorem  is  proved. 

n—*oo 

Remark  In  practice  it  is  frequently  the  case  that  the  admissible  parameter  set  Q  is 
also  infinite  dimensional;  for  example,  when  the  unknown  parameters  to  be  identified  are 
elements  in  a  function  space.  In  this  situation  the  set  Q  must  also  be  discretized  and 
Theorem  4.1  must,  to  a  certain  degree,  be  modified.  For  each  m  =  1,2, . . .  let  Im  :  Q  C 
Q  -+  Q  be  a  continuous  map  with  range  Qm  =  Im(Q)  in  a  finite  dimensional  subspace 
of  Q  and  with  the  property  that  lim  Im(q )  =  9,  uniformly  in  9  for  9  e  Q.  We  then 
define  the  doubly  indexed  sequence  of  approximating  identification  problems  {ID™)  by 
letting  problem  {ID™)  be  the  problem  {IDn)  with  the  set  Q  replaced  by  the  set  Q™.  The 
modification  to  Theorem  4.1  would  state  that  each  of  these  problems  admits  a  solution 
9J?  G  Q™  and  that  the  sequence  {9^}  will  have  a  Q-convergent  subsequence  whose  limit 
is  in  Q  and  is  a  solution  to  problem  {ID).  We  note  that  each  of  the  problems  {ID™) 
involves  the  minimization  of  a  continuous  functional  over  a  compact  subset  of  Euclidean 
space  subject  to  finite  dimensional  state  space  (ODE)  constraints.  As  such,  for  each  m  and 
n,  problem  {ID™)  may  be  solved  using  standard  computational  algorithms  and  techniques. 
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JL  An  Exmni?l<? 


In  order  to  illustrate  the  application  of  the  general  theory  that  was  presented  above,  we 
consider  the  problem  of  estimating  or  identifying  the  nonlinear  damping  term  in  the  forced 
one  dimensional  wave  equation 

(M)  (M)  +  l>(9)  (<»»/))  -  (M)  9  /  (M)  .<  >  0,  0  <  77  <  1 

with  the  Dirichlet  boundary  conditions 

(5.2)  u(f,0)  =  0,  u(f,l)  =  0,  t>0 

and  initial  data 

(5.3)  w(0 ,77)  =  u°(rj),  -^-(0,7)  =  u\v),  0  <  77  <  1. 

We  note  that  for  definiteness  we  have  chosen  the  Dirichlet  boundary  conditions  (5.2),  how¬ 
ever  all  of  the  discussion  below  would  remain  with  any  of  the  usual  self-adjoint  boundary 
conditions  (i.e.  Neumann,  Robin,  etc.).  We  assume  that  a  >  0,/  £  Lz((0,T)  x  (0,  l)),u°  £ 
Hq{ 0, 1)  and  u1  £  £2(0, 1).  Let  Q  be  a  metric  space  and  let  Q  C  Q.  be  compact.  For  each 
q  £  Q  we  assume  that  the  mapping  b(q)(-)  :  R  — +  2*  satisfies  the  following  conditions: 

(bl)  The  set  {<p  £  H]  (0,1):  there  exists  a  t/>  £  H~x  (0,1)  such  that  tp(rj)  £  b(q)(<p(r])),  a.e. 
tj  £  (0, 1)}  is  independent  of  q  £  Q,  and  0  £  b(q)( 0), 

(b2)  The  mapping  q  — »  b(q){9)  is  lower  semi-continuous  from  Q  C  Q  into  2®  for  almost  every 

6£R, 

(b3)  For  each  q  £  Q  the  mapping  b(q)( •)  is  nondecreasing  and  for  some  A  >  0  the  inclusion 
6  -f  \b(q)(9)  3  £  has  a  solution  9  £  R  for  each  £  £  R  (i.e.  b(q)(-)  is  maximal  monotone 
in  R), 

(b4)  There  exists  a  polynomial  p,  independent  of  q  £  Q  for  which  |0|  <  p(\9\)  for  all  9  £ 
6(<j) (|0|),  and  almost  every  9  €  R. 
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To  reformulate  the  initial-boundary  value  problem  (5.1)  ••  (5.3)  as  an  abstract  system  of 

the  form  (3.1),  (3.2),  we  set  H  —  £2(0,1)  endowed  with  the  standard  inner  product,  and 

1  2 

let  V  =  Hi  (0,1)  be  endowed  with  the  norm  ||</j||  =  f  \Dy>{yi)y  drj.  In  this  case  we  have 

0 

V*  =  H~l  (0,1)  together  with  the  dense  and  continuous  embeddings  V  <— *  H  <— ►  V*.  We 
define  the  operator  A  €  C(V,  V *)  by 

1 

(5.4)  <A<p,ip  >=  a  J  Dtp  (ri)  Dif>  (ri)  dri,  e  Hq  (0,1) . 

0 

It  is  easily  shown  and  well  known  that  the  operator  A  given  by  (5.4)  satisfies  conditions 
(Al)  -  (A4).  (Note  that  in  the  present  example  the  operator  A  is  assumed  known  and  as 
such  does  not  depend  on  q.) 

To  define  the  operator  B(q )  we  follow  the  treatment  in  [18]  and  use  the  notion  of  a 
subdifferential  of  a  proper  convex  lower  semicontinuous  mapping  (see  [10]).  For  each  q 
6  Q  let  bo(q)(-)  denote  the  minimal  section  of  the  mapping  b(q)(-).  That  is  b0(q)(-)  is  the 
single-valued  mapping  from  R  into  R  defined  by  bo(q)  =  9,  where  9  is  the  unique  element 
in  b(q)(9)  of  minimal  absolute  value.  That  the  minimal  section  of  b(q)( •)  is  well  defined  is 
a  consequence  of  condition  (b3).  Since  Dom(6( <?)(•))  =  R,  the  proper,  convex,  lower  semi 
-continuous  function  j(-;q ) :  R  — »  R  can  be  defined  by 

9 

3  (0;$)  =  J  bo(q)(Odt 

0 

where  R  denotes  the  extended  real  numbers.  It  can  be  shown  that  j(-;  q)  is  bounded  below 
by  an  affine  function  so  that  j(9;  q)  >  -00,  9  £  R,  and  that  dj  (•;  q)  —  b(q)  (•),  where  d 

denotes  the  subdifferential  operator.  For  each  q  G  Q  define  7  (-;g)  :  Hq  (0, 1)  — ►  R  by 

7(^4)  =  /  if  ;(v>  (•);?)  €  £1  (0,1) 

(  +00  otherwise. 

Then  7 (*;g)  is  also  proper,  convex  and  lower  semi-continuous  and  we  define  B(q)  : 
Dom(5)  CV  ->2V‘  by 

(5.5)  B  (q)  v?  =  dy  (ip;  q) , 
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where  <p  G  Dom(B(7))  =  {p  G  Hq  (0, 1) :  d~f(tp\ q)  ^  0}.  (The  fact  that  Dom(5(f/))  = 
Dom(B)  is  independent  of  q  is  a  consequence  of  condition  (bl).  )  It  can  be  shown  (see,  for 
example,  [10,  p.  61])  that  for  <p  G  Dom(B),  rp  G  B{q)<p  if  and  only  if  rp(rj)  G  b(q)(<p(r))), 
a.e.  r\  €  (0, 1),  and  via  conditions  (bl)  -(b4),  that  the  operators  B{q)  given  by  (5.5)  satisfy 
conditions  (Bl)  -  (B4). 

To  illustrate  the  formulation  of  an  inverse  problem  we  take  the  observation  space  Z  to 

If 

be  given  by  Z  —  x  {R*  x  Li  (0, 1)}  and  the  performance  index  from  C(0,T;  V  x  H)  into 
the  nonnegative  reals  to  be  given  by  the  weighted  least-squares  functional 

V  t 

«=1  j=l 

for  (u,u)  G  C(0,T;V  x  H)  and  z  =  ((2}, z?) , .  •  •  G  Z  with  pi,&i  >  0 ,i  = 

1,2,...  , v  ,  0  <  <1  <ti<  ...  <tv  <T,  and  0  <  rji  <  772  . . .  <  qt  <  1. 


J  \u{Uyr))-z]  (t7)|2cM 


We  provide  two  specific  examples  of  possible  choices  of  the  parameter  space  Q,  the 
admissible  parameter  set  Q ,  and  mappings  &($)(■)  :  R  — ►  2R  satifying  conditions  (bl)  - 
(b4).  A  relatively  simple  example  might  involve  the  estimation  of  the  const  ant.  parameters 
q  =  ( oto,Po,00 )  in  the  saturating  polynomial  function 


Km  = 


0o  |^|a°  sgn{6) , 
0o  |flo|tto  sgn(9), 


-60  <0<  d0 

\e\  >  e0. 


In  this  case  we  would  have  Q  =  R3  and  Q  =  {(ao,/?O)0o)  :  0  <  o0  <  d,0  <  0o  <  0,0  < 
do  <  6}  for  some  d,  0, 6  >  0  given  and  fixed.  We  note  that  when  <*0  =  0,  in  order  to  satisfy 
the  maximality  in  condition  (b3)  we  must  take 


We  note  also  that  when  d  <  2  we  have  Dom(B)  =  V  and  that  the  operators  B(q)  given 
by  (5.5)  map  Jlf-bounded  subsets  of  V  into  V*-bounded  subsets,  uniformly  in  q  for  q  G  Q. 


Thus  the  hypotheses  of  Theorem  3.2  are  satisfied  and  we  see  that  Dom (A(q))  =  VxH. 
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A  second,  somewhat  more  challenging  problem  might  involve  the  estimation  of  a  func¬ 
tional  parameter.  In  this  case  we  take  C  s=  Cfl(R),  the  set  of  all  bounded,  continuous,  real 
valued  functions  on  R  endowed  with  the  supremum  norm.  We  take  Q  to  be, the  (2-closure 
of  the  set 

{q  G  Q  :  q  (6)  =  -q(-B),6q(0)  >  0,for  0  G  R, 

(5.6)  |?(0)|  =9(^0), for |0|  >  90}q  G  H 1  (-0oA)» 

i9i«‘(-flo,«o)  ^  K°'  v'w  -  °* for  a-c- 6  G 

where  0o  and  Ko  are  given  positive  constants.  It  is  not  difficult  to  argue  (via  the  Arzela- 
Ascoli  Theorem)  that  Q  is  a  compact  subset  of  Q.  For  q  G  Q  we  set  b(q){6)  =  q(0),0  G  R. 


With  regard  to  illustrating  the  approximation  or  discretization  of  the  state,  we  briefly 
outline  a  linear  spline-based  scheme.  For  each  n  =  1,2, . . . ,  let  Hn  =  span{y>3n}  j-\  where 
is  the  j  -  th  linear  B-spline  on  the  interval  [0, 1]  defined  with  respect  to  the  uniform 


mesh  {0, 1/n,  2/n, . . .  ,1).  That  is, 

Vn  (V)  =  ' 


0 

nr)  -  j  +  1 
;  +  1  -  nrj 
0 


ill  <  n  <  i- 

n  —  /  —  n 

i  <  „  <  i±i 

n  —  '  —  n 
^  <«?<!, 


;  =  1,2, .. .  ,n  -  1.  It  is  immediately  clear  that  Hn  C  V  =  Hq(0,  1),  for  each  n.  Let 
Pn  :  H  — +  Hn  denote  the  orthogonal  projection  of  £2(0, 1)  onto  Hn  computed  with  respect 
to  the  usual  L2  inner  product.  Using  well  known  estimates  for  interpolatory  splines  (for 
example,  those  found  in  [23])  it  can  be  argued  that  condition  ( P )  is  satisfied. 


In  the  example  involving  the  identification  of  a  functional  parameter,  the  admissible 
parameter  set  Q  chosen  to  be  the  Q-closure  of  the  set  given  in  (5.6)  can  be  discretized  as 
follows.  For  each  m  =  1,2,...  and  0  G  R+  let  {tpj1  (•, @)} JL0  denote  the  standard  linear 
J3-splines  on  the  interval  [0,0]  defined  with  respect  to  the  uniform  mesh  {0,  . . .  ,  0} 

and  then  extended  to  a  continuous  function  on  the  entire  positive  real  line  via  (0;  0)  = 
rPJ1  (0;  0)  ,0  >  0.  For  q  G  Q  set 

m 

(5.7)  (/”■«) (*)  =  £«(&)  *T  (l9l i #.) " m 

>-l 
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for  9  €  R  whore  9q  is  the  number  in  (O,0o)  for  which  |<7(0)|  =  q{9q) ,  |#|  >  <V  (Note  that 
the  lower  limit  in  the  sum  in  (5.7)  is  1  rather  than  0  since  q  G  Q  implies  q{ 0)  =  0.)  Using 
the  Pcano  Kernel  Theorem  (see  [21]  p.  22)  it  can  be  argued  that 

-  ?lloe  =  S«P  I ~  ?!  <  \  (S)’/2  A'o 

and  consequently  that  lim  Imq  =  o,  uniformly  in  q  for  q  €  Q .  We  would  then  set 

m— *oo 

Qm  =  Im(Q). 


For  Hn  and  Qm  as  defined  above,  the  finite  dimensional  initial  value  problem  (4.1),  (4.2) 
takes  the  form 

(5.8)  Mnwn{t)  +  Cn(wn(t)]qm)  +  Knwn(t)  =  Fn(t),  0  <  t  <  T, 

(5.9)  Mnwn( 0)  =  w°n,  Mn u>„(0)  =  wl, 

where  F„(<),u>°,  and  w\  are  the  (n-l)-vectors  whose  i-th  components  are  given  by  F'n(t) 
—  <  /(Oi^n  >iwn  =<  >  and  =  <  u1,^  >,  respectively,  and  M„  and  Kn 

are  the  (n  —  1)  x  (n  —  1)  matrices  whose  (ij)-th  entries  are  given  by  M =<  ip'n,<p]n  > 
and  =<  a<p*n,(pJn  >.  The  vector  function  Cn( •  ;  qm )  :  Rn_1  — ►  Rn_1  is  given  by 

r? 

Ch(v,qm)=  J  {nx  —  i -f  1}  ({nx  —  i}  {t>*  —  v*-1 }  +  u‘)  dx 

i-i 

n 

n 

+  J  {i  +  1  -  nx)  qm  ({nx  -  i}  {u,+1  —  u*}  -f  v')  dx, 

1 

n 

i  —  1,2,  •  •  •  ,n  —  1,  for  v  €  Rn_1  with  v°,vn  =  0.  If  wn{ •  ;gm)  is  the  solution  to 

the  second  order  initial  value  problem  (5.8),  (5.9)  corresponding  to  qm  €  Qm,  then 

«n(U9m)=  E  <(<;?m)^  andtin(<;?m)  =  f)  w}n(t;  qm)v3n,  for  t  €  [0,T].If9m  6  Qm 
>=i  j= i 

m 

is  given  by  qm{9)  =  £  AjP^jPO^Ij^V")  S0n(0),  ^  €  R,  the  identification  problem  {ID™) 

>*1 
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becomes  one  of  determining  parameters  (g™,  i®qm)  >n  some  compact  subset  of  Rm+1 
which  minimize  $n(?m)  =  <£(un(gm),  un(qm)\  z). 


A  discussion  of  some  implementation  questions  relevant  to  the  schemes  outlined  (in 
particular  with  regard  to  supercomputing)  together  with  numerical  results  can  be  found 
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